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O Abstract 

A /3-skeleton, /3 > 1, is a planar proximity undirected graph of an Euclidean point 

set where nodes are connected by an edge if their lune-based neighborhood contains 

I— I no other points of the given set. Parameter (3 determines size and shape of the 

W nodes' neighborhoods. In an excitable /3-skeleton every node takes three states — 

W resting, excited and refractory, and updates its state in discrete time depending on 

P3 states of its neighbors. We design families of /3-skeletons with absolute and relative 

J-H thresholds of excitability and demonstrate that several distinct classes of space-time 

. ^. excitation dynamics can be selected using /3. The classes include spiral and target 

waves of excitation, branching domains of excitation and oscillating localizations. 

^ Keywords: proximity graphs, /3-skeletons, excitation, waves, localizations, space-time 

"^ dynamics, pattern formation 

o 

o 

^ 1 Introduction 



^ 



A planar graph consists of nodes which are points of Euclidean plane and 
edges which are straight segments connecting the points. A planar proximity 
^ graph is a planar graph where two points are connected by an edge if they 

are close in some sense. Usually a pair of points is assigned certain neighbor- 
hood, and points of the pair are connected by an edge if their neighborhood 
is empty. Delaunay triangulation [8j, relative neighborhood graph [10] and 
Gabriel graph [tl6j, and indeed spanning tree, are most known examples of 
proximity graphs. /5-skeletons, proposed in fT2j, form a unique family of prox- 
imity graphs monotonously parameterised by parameter /3. 

Proximity graphs found their applications in fields of science and engineerings: 
geographical variational analysis ^9lll6ll20j . evolutionary biology [15], spatial 
analysis in biology [TSllGllTlfTT] , simulation of epidemics [23j . Proximity graphs 
are used in physics to study percolation [5] and analysis of magnetic field [22] . 
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Engineering applications of proximity graphs are in message routing in ad hoc 
wireless networks, see e.g. p!4ll2 111911171124) . and visualisation [18j. Road net- 
work analysis is yet another field where proximity graphs are invaluable. Road 
networks are well matched by relative neighborhbood graphs, see e.g. study 
of Tsukuba central district [251126] . Biological transport networks also bear re- 
markable similarity to certain proximity graphs. Foraging trails of ants [1] and 
protoplasmic networks of slime mold Physarum polycephalum [S^'S] are most 
striking examples. 

Structure of proximity graphs represents so wide range of natural systems 
that it is important to uncover basic mechanism of activity propagation on 
the graphs, which could be applied in future studies of particular natural 
systems. This is why we undertook computational experiments with excitable 
/5-skeletons to check how space-time dynamics of excitation on /3-skeletons 
depends on /3. The paper is structured as follows. In Sect. [2] we introduce 
models of excitable /3-skeletons. Phenomenology of skeletons with absolute 
threshold of excitation (a node excites depending on an absolute number of 
its excited neighbors) is provided in Sect. [3J Space-time dynamics of skeletons 
with relative threshold of excitation (a node excites depending on a ratio of 
excited neighbors) is studied in Sect. [4| Results of computational experiments 
are summarized in Section [H 



2 The model 



Given a set V of planar points, for any two points p and q we define /3- 
neighborhood Uf3{p^q) as an intersection of two discs with radius P\p — q\/2 
centered at points ((1 — 1)]9, |g) and (|]9, (1 — |)g), (5 > 1 [12ul0j . see examples 
of the lunes in Fig. [T] Points p and q are connected by an edge in ;5-skeleton 
if the pair's /3-neighborhood contains no other points from V. 

A /3-skeleton is a graph G^iy) = (V,E,/5), where nodes V C R^, edges E, 
and for p, g G V edge (pq) G E if [/^(p, q) H V/jp, q} = 0. Parameterization 
/5 is monotonous: if /5i > /32 then G^^(V) C G/^sfv) [T2l[T0] . A /5-skeleton is 
Gabriel graph [16j for ;5 = 1 and the skeleton is relative neighbourhood graph 
for ;5 = 2. We consider only skeletons with 1 < /5 < 2 because ;5-skeletons are 
non-planar for ;5 < 1 and they are disconnected for /5 > 2. 

We study ;5-skeletons of planar set of n discs. Centers of the discs form set 
V. Each disc has a radius 2.5 units. The discs are randomly distributed in 
a large disc with radius 480 units. We undertake computational experiments 
with n varying from 1000 to 15000 (Fig. [2]). Number of nodes per se is not as 
important as density of their packing therefore we refer to any particular V 
via density of node packing ip = 0.027, • • • , 0.407. 
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Fig. 1. Examples of the lunes between two points (small circles) for various values 
of/3. 



An excitable /3-skeleton is defined as follows. A node p G V is a finite state 
machine. Every node takes three states: resting (o), excited (+) and refractory 
(— ). A node updates its states in discrete time depending on states of its 
neighbors. All nodes update their states simultaneously. 

We assume that a resting node excites depending on a number of excited 
neighbors. If a node is excited at time t the node takes refractory state at 
time step t + 1, independently on states of its neighbors. Transition from 
refractory to resting state is also unconditional. 

Let iy{p) = {g G V : (pq) G E} be a neighborhood, or set of neighbors, of node 
]9, p^ a state of node p at time step t, (J^{p) a number of excited neighbors of 
p at step t, and d{p) a degree, or a number of neighbours |^(p)|, of node p. 
Then the node-state transition function can be defined as follows: 



P' 



t+i 



f+, ifCiuipY) 

-, iip^ = + 

o, if p^ = — 



(1) 



We consider two versions of the excitation condition C{v{py) G { True, 
False }: 





(a) /3 = 1 



(b) /3 = 1.5 




(c) /3 = 1 



(d) (3 = 1.5 



Fig. 2. Examples of /3-skeletons: (a) n = 100, (f = 0.0027, and /3 = 1, (b) n = 100, 
(^ = 0.0027, and /3 = 2, (c) n = 5000, (p = 0.136, and /3 = 1.5, (d) n = 10000, 
(^ = 0.271, and/3 = 1.5. 

• Absolute excitability: C{v{py) = cr^(p) > ^, ^ = 1, • • • 

• Relative excitability: C{u{py = (^7^^ > e), e G [0, 1] . 

Absolute excitability is the most common approach of defining rules in ex- 
citable discrete systems however it does not account for diversity of node 
degrees in disordered systems. Thus we also explore skeletons with relative 
excitability. 



In the paper we illustrate space-dynamics of excitation by snapshots of skele- 
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Fig. 3. Activity a versus /3 for various densities (f of node packing in /3-skeletons with 
absolute excitation threshold = 1. Skeletons are developed from initial random 
configuration, where a node is assigned excited state probability 0.1 and is resting 
otherwise. 

ton configurations. We analyze integral dynamics of skeletons using activity a. 
The activity a is a ratio of excited nodes to a total number of nodes averaged 
fixed number time steps. The activity a is measured after initial transient pe- 
riod, when excitation patterns are given a chance to occupy the whole skeleton 
(usually a hundred time steps is sufficient). 



3 Absolutely excitable skeletons 



In /5-skeletons governed by absolute excitability rule excitation persists only 
for threshold 9 = 1. 

Finding 1 Activity a is proportional to density Lp and inversely proportional 
to f3. 

Figure [3] shows exactly how activity a depends on /3. Note that maximum 
possible a equals 0.33, for sustainable activity, because an excited node always 
becomes refractory, and a refractory node always becomes resting. Thus even 
if the whole skeleton is active, only third of its nodes are excited at any given 
time step. The activity decreases polynomially with increase of /5, degree of 
the polynomial is proportional to density (f of nodes in the skeleton. 

Spatial counterparts of integral activity are shown in Fig. [4} These are exam- 
ples of excitation patterns emerged after a single node of the resting skeleton 
was excited. In ;5-skeletons with high density of nodes packing Lp = 0.271 and 
0.407 excitation of a single node of a resting skeleton leads to formation of 
a generator of spiral and target waves if 1 < ^ < 1.4. The most noticeable 
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Fig. 4. Configurations of excitation on ^-skeletons with node density 
(^ = 0.027, 0.136, 0.271, 0.407 and absolute threshold of excitation 6> = 1. All config- 
urations are developed from a localized initial excitation: all nodes but are resting, 
one node is excited. The configurations are snapped after 50 steps of development. 
Red (light-gray) coloured node are excited, blue (dark-gray) coloured nodes are 
refractory. 




Fig. 5. Dependence of critical values jid and j3o on ^p. Data points for jid are shown 
by circles, for ^o by triangles. Trend lines are ^^ — l.Q^Lp + 1.05, R^ — 0.961 and 
Pa = l-42(/p + 1.45, R^ = 0.964. 



pattern of the waves is observed for Lp = 0.407 and /3 = 1. Increase of /3 leads 
to break up wave- fronts and only localized compact clusters of activity remain 
on the skeleton. 

What values of (p and /3 cause excitation fail to span the whole graph? For what 
values of (p and /5 only tiny localized clusters of oscillating activity remain in 
a skeleton? We call values /3 critical if patterns of excitation change abruptly, 
e.g. excitation dynamics changes from target waves filling the whole skeleton 
to slowly growing, /3d^ or even localized domains of activity, oscillators, /3o. 

Finding 2 Critical values /3d and /3o are linearly proportional to density ip of 
nodes in skeleton. 

Assuming that a skeleton fully occupied by target waves has activity level 
a = 0.33 even tiniest decrease in a indicates presence of a stationary resting 
domain. See example in Fig. [4| 99 = 0.271 and /5 = 1.4. Let us assume that 
excitation fails to span a skeleton when a < 0.32 and only bounded domains 
of excitation are present if a < 0.25. Let /3d and /3o be critical values of /3 thus 
that if for a given /5-skeleton value /3 exceeds /3d then the skeleton is not fully 
occupied by excitation, and if /3 > /3o then only localized oscillators are formed. 
Figure [5] presents critical values of /3 computed for 25 sample configurations 
of /3-skeletons with Lp varying from 0.027 to 0.407. What are structural corre- 
lates responsible for the transition from a wide-spread excitation to localized 
domains? 



Finding 3 Excitation ceases to occupy a whole skeleton when the mode of 
node degree distribution of the skeleton changes from 4 to 3. 
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Fig. 6. Mode (a) and average (b) node degree of /3-skeletons for /3 G [1,2] and 
(f = 0.027, • • • , 0.407. In (a) entries are boldfaced to increase visibility. In (b) entries 
corresponding to /3o((f) are underlined. 
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Fig. 7. Parameterisaton of a space-time dynamic of /3-skeleton with absolute exci- 
tation threshold by node density ^p and j3. 



We found this by directly comparing /3^ and integral characteristics of node 
degree distributions for a range of /3 and ip (Fig. [6]). Mode of node degree 
distribution is not helpful however in detecting skeletons supporting localized 
oscillators. As shown in Fig. [6] skeletons with mode 3 of degree distribution can 
exhibit very large and very small domains of oscillating activity. An average 
node degree of a skeleton gives us a bit more detailed picture of 99-/3-induced 
structural transitions and thus can be used to detect when localized oscillators 
are formed. 



Finding 4 A skeleton exhibits only localized oscillating domains of activity 
when average node degree drops below 3. 

The average node degrees corresponding to onset of localized oscillations are 
underlined in Fig.JGJ Further comparison of average node degrees (Fig.JGJ), and 
corresponding space-time configurations of activity (see e.g. Fig. ^ allows 
us to state the following observation. Average node degree 5 of a skeleton 
determines space-time dynamics of excitation as follows: 

• 3.9 < 8\ classical wave-like dynamics, generators of target waves, ordered 
dynamics (Fig. JTJA), 

• 3.7<5<3.9: generators of wave- fragments, disordered excitation activity 

(Fig.|7^), 

• 3.0 < 5 < 3.7: formation of slowly propagating domains of excitation, some 



loci of a skeleton remain resting (Fig. [7|C), 
• 5 < 3.0: only small localized stationary domains of oscillating activity are 
formed (Fig. [7|D). 



4 Relatively excitable skeletons 



Here we consider behavior of skeletons with node density Lp = 0.407. Skele- 
tons with density (p = 0.407 (15,000 disc-nodes) show a wide representation 
of structural properties because they exhibit widest range of modes of node 
degree distributions (Fig. [6]). In computational experiments we perturbed a 
resting skeleton with small excitation: every node is assigned an excited state 
with probability 0.1. We allowed the skeleton to develop its excitation patterns 
for 100 iterations and recorded activity a then. Examples of integral activity 
of skeletons for /3 = 1, 1.5 and 2 are shown in Fig. [8| We see that drastic 
changes in activity levels occur when e changes from 0.199 to 2 (first drop in 
activity level is most clear in Fig. |8^), and when e changes from 0.249 to 0.25 
(second drop in activity level, see e.g. Fig. (s): 



More detailed description of how a depends on /5 is provided in Fig. |9] For 
e G [0,0. 2[ activity level a is not changed with increase of /3, 1 < /3 < 1.7. 
However as soon as /3 reaches 1.8 activity level slightly decreases. The activity 
level drops down by one third when /3 = 2 (Fig. [9]). For e G [0.25, 0.33[ activity 
monotonously increases with increase of /3. For e G [0.2, 0.249[ activity a grows 
with growth of /3 till /5 = 1.6. The activity remains unchanged for ;5 = 1.7 and 
then decreases with further increase of /3. Let us discuss space-time dynamics 
of these /3-skeletons for selected values of /3. 



4.1 /3 = l 



For e G [0,0.166] skeletons exhibit classical properties of discrete excitable 
media. A single excitation causes formation of a spiral wave or a generator 
of spiral or target waves. Successions of circular waves fill the whole skeleton 



(Fig. lO^ib). As soon as e exceeds 0.166 excitation wave- fronts break up into 



separate wave fragments. Thus branching domains of excitation activity are 



formed (Fig. 10:). 



Only tiny oscillating activity domains emerge for e > 0.2 (Fig.jlOp). Life-cycles 



of three oscillators are shown in Fig. 11, Left oscillator in Fig. ^Tlis a typical 



one. Excitation wave, consisting of one excited state and one refractory state, 
runs along a cycle of five nodes. The oscillator thus has period five (Fig.[TT^-f). 
The oscillator conserves number of non-resting states. 
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(b) /3 = 1.5 



0.4 



0.3 



0.2 

a 

0.1 



0.0 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 

8 

(c) /3 = 2 

Fig. 8. Activity a versus relative excitation threshold e of /3-skeletons for (a) P — 1, 
(b) P = 1.5, (c) /5 = 2. 
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Fig. 9. Dependence of activity level a of /3-skeletons on /3 and e. Initially resting 
skeleton is set to a random configuration of excitation by exciting each node with 
probability 0.1. Values of a are measured 100 iterations after initial excitation of 
skeletons. 



Oscillator in the top part of snapshots in Fig. [TTJhas period three. It changes 
its state between configurations with two excited and one refractory and one 
excited and two refractory states. The activation does not cycle but is ex- 
changed between several nodes. First one node a excites its two neighbors. 
These excited nodes convey excitation to two other nodes (their neighbors), 
which in turn excite the original node a (Fig. [TT^-d). 



Oscillator shown in right part of snapshots in Fig. [TT] consists of four excited 
and four refractory states, the number of the non-resting states is conserved. 
The oscillator has period three (Fi g. [TT^ -d). This oscillator behaves like a 
breather, with its compressed (Fig. [lip), intermediate (Fig. [TI^) and fully 
expanded (Fig. [TT|3) configuration. 



Skeleton becomes non-excitable, i.e. no excitation persists, when e exceeds 
0.25. 
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(c) /3 = 1, 6 = 0.167 







(b) /3 = 1, e = 0.15 





(d) /3 = 1, e = 0.2 




(e) /3 = 1, e = 0.24 

Fig. 10. Snapshots of excitation dynamics in /3-skeletons (/3 = 1) with relative 
excitation threshold, (a)-(d) singular initial excitation, (e) random initial excitation, 
a node is excited with probability 0.1. Red (light-gray) coloured node are excited, 
blue (dark-gray) coloured nodes are refractory. 
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Fig. 11. Examples of oscillators in /3-skeletons (/3 = 1) with relative excitation 
threshold e = 0.24. Red (light-gray) coloured node are excited, blue (dark-gray) 
coloured nodes are refractory. 



12 (3 = 1.5 



Propagating domains of excitation are circularly shaped however they do not 
show any pronounced wave-fronts, the excitation patterns are rather disor- 
dered for e G [0,0.199] (Fig. [T2^b). With e exceeding 0.199 activity domains 
change their shapes form circular to branching, tree-like propagating domains 



(Fig. 12:). Size of such branching domains decreases with increase of e. When 



e > 0.25 all domains cease propagating and just tiny clusters, oscillators, of 



activity are formed (Fig. 121). A typical oscillator is a cycle of three or four 
nodes around which a quasi-one-dimensional excitation propagates. Sometimes 
the running excitation sends 'sparks' of activity to lateral nodes; these sparks 
extinguish in few steps of development. No excitation persists in skeletons for 
e > 0.33. 
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Fig. 12. Snapshots of excitation dynamics in /3-skeletons (/3 = 1.5) with relative 
excitation threshold. Red (light-gray) coloured node are excited, blue (dark-gray) 
coloured nodes are refractory. 

13 /5 = 2 



Single excitation gives birth to slowly propagating irregularly-shaped domains 



of excitation activity (Fig. 13 a.) . The domains originated form a single excita- 
tion never span the whole skeleton. Sizes of domains decrease with increase of 
e(Fig.[l3bc). 



Only tiny clusters of excitation are formed for e G [0.25,0.333[ (Fig. [l3|d). 
A life-cycle of most typical oscillator is shown in Fig. [14] This oscillator has 
period seven. Its core structure is a singleton-excitation (accompanied by a 
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Fig. 13. Snapshots of excitation dynamics in /3-skeletons (/3 = 2) with relative 
excitation threshold. Red (light-gray) coloured node are excited, blue (dark-gray) 
coloured nodes are refractory. 

refractory tail) running around seven-node cycle, anti-clockwise. At certain 
moments of its life-cycle the excitation excites few nodes neighboring to the 



'core cycle'. Thus sparks of excitation are formed (Fig. 14:^f). 

No excitation persists for e > 0.333. 

Finding 5 Excitable beta-skeletons with relative threshold of excitation ex- 
hibit the following classes of space-time dynamics: (1) spiral and target waves 
similar to classical excitation waves, (2) omnidirectional propagating patterns 
of disordered excitation, (3) localised domains of excitations, (4) branching 
domains, and, (5) tiny oscillators. 
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(a)t 



(b) t + 1 



(c) t + 2 



(d) t + 3 




(e) f + 4 



(f ) i + 5 



(g) t + 6 



(h) f + 7 



Fig. 14. Examples of oscillators in /3-skeletons (/3 = 2) with relative excitation 
threshold e = 0.29. Red (light-gray) coloured node are excited, blue (dark-gray) 
coloured nodes are refractory. 
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Fig. 15. Parameterisaton of space-time dynamics of /3-skeletons with relative exci- 
tation threshold by /3 and e. 



Position of the classes on f3 — e-plane is shown in Fig. [T5j There are sharp 
boundaries between classes, characterized by abrupt changes in morphology 
of excitation patterns. Critical values of relative excitability threshold are 
e = 0.167, 0.2, 0.25 and 0.33. These thresholds reflect situations when a resting 
node excites if at least one of six (e = 0.166), one of flve (e = 0.2), one of 
four (e = 0.25) and one of three (e = 0.33) neighbors is excited. Excitability 
threshold e = 0.166 is critical only for ;5-skeletons with /3 very close to 1. This 
is because such skeletons have mode 5 of node degree distribution and average 
node degree 5.1 (Fig. |6| 99 = 0.407), i.e. nodes with six neighbors are second 
(after nodes with flve neighbors) dominating nodes on the graph. For similar 
reasons e = 0.2 is not critical threshold value for /3-skeletons with /5 close to 
2 because over half of the nodes in such skeletons have three neighbors, and 
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(a) 



0.204 




(b) 



0.199 






(c) e = 0.204 



(d) e = 0.25 







(e) e = 0.249 



(f) e = 0.25 



Fig. 16. Stability of oscillators under changes of e: (a)-(c) /3 = 1, (d)-(f) /3 = 2. 



nodes with two or four neighbors are second dominants in these skeletons. 



4^4 Stability of localized oscillators 



Finding 6 Localized excitations^ or oscillators^ are stable under changes of 
relative excitability threshold. 



The finding is illustrated for /3 = 1 and /3 = 2 in Fig. [T6| We choose values 
of threshold e at the edge of excitability — e = 0.204 for /3 = 1 and e = 
0.25 for /3 = 2 — and perturb /3-skeleton with a random configuration of 
excitations. Few localized oscillators emerge (Fig. [iG^d). Then we decrease 
excitability thresholds — dynamically, without resetting the skeleton to resting 
state ^ to 6 = 0.199 {/3 = 1) and e = 0.249 (/5 = 2). The oscillators start 
to generate excitation wave- fragments. The wave-fragments emitted by the 
oscillators fill the skeletons (Fig.[T6)3e). When we reduce excitability thresholds 
back to their original levels 0.204 and 0.25 traveling wave- fragments disappear 
but oscillators remain untouched, they stay where they were before relative 



excitability threshold changed (Fig. 16 :f). 
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5 Discussion 



We studied discrete excitation dynamics on a smoothly parameterised family 
of planar proximity graphs — /3-skeletons. The skeletons were constructed on 
a set of uniform discs packed into a larger disc. Thus each skeleton is char- 
acterised by /3 and density of node packing. We considered rules of absolute 
excitability (a resting node excites if number of its neighbors exceeds certain 
threshold) and relative excitability (a resting node excites if a ratio of its 
excited neighbors to a total number of neighbors exceeds certain threshold) . 

In computational experiments we found that overall level of activity in an 
absolutely excitable /5-skeleton is proportional to node packing density and 
inversely proportional to /5. We demonstrated that space-time dynamics of 
absolutely excitable /5-skeletons can be classified as follows: spiral and target 
waves, disordered dynamics with irregularly traveling wave- fragments, slowly 
spreading domains of excitation, and small localized stationary domains of 
oscillatory activity. Transition of a skeleton between these classes is controlled 
by changing /3 and density of node packing. Both /3 and density affect distri- 
bution of node degrees. We provided a parameterization of the classes using 
just one parameter — average node degree. 

There are five classes of space-time activity of the relatively excitable /3- 
skeletons with tightly packed nodes (density 0.407, 15000 nodes): spiral and 
target waves, omnidirectional propagating patterns of disordered excitation, 
localized domains of excitations, branching domains of activity and tiny oscil- 
lators. The classes are controlled by /3 and the relative excitability threshold 
e. Tiny oscillators usually emerge just below edge of excitability. The tiny 
oscillators are stable under relative excitability threshold. Decreasing thresh- 
old we can increase number of excited loci around the oscillators but we can 
not modify oscillators themselves. As soon as we increase the threshold all 
excitation activity but original oscillators disappear. 

The class of branching domains is a transitional class occupying a part of 
/3 — 6-space between full excitability and non-excitability. With regards to 
propagating patterns, waves are pronounced for /5 = 1, and are alike classical 
excitation waves in discrete media. Increase of /5 to 1.5 causes wave- fronts to 
break up into separate wave- fragments. 

We believe our results will find their applications in analysis of spreading 
crowd activities, e.g. riots, on a city's streets; distributed containment of traffic 
jams on motorways networks; study of foraging patterns of social insects and 
myxomycetes; control of excitation dynamics and propagation of defects in soft 
matter; design and simulation of conglomerates of simple excitable elements. 
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